Process reconstruction from incomplete and/or inconsistent data 
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We analyze how an action of a qubit channel (map) can be estimated from the measured data that 
are incomplete or even inconsistent. That is, we consider situations when measurement statistics is 
insufficient to determine consistent probability distributions. As a consequence either the estimation 
(reconstruction) of the channel completely fails or it results in an unphysical channel (i.e., the 
corresponding map is not completely positive). We present a regularization procedure that allows 
us to derive physically reasonable estimates (approximations) of quantum channels. We illustrate 
our procedure on specific examples and we show that the procedure can be also used for a derivation 
of optimal approximations of operations that are forbidden by the laws of quantum mechanics (e.g., 
the universal NOT gate). 



I. INTRODUCTION 

For any reliable quantum information processing it 
is important to know how states of quantum systems 
are transformed under the action of quantum channels 
(maps). It is therefore essential to develop tools by means 
of which we can acquire knowledge about properties of 
quantum channels. Providing we have no prior knowl- 
edge about an action of a particular transformation our 
task is to determine characteristics of a corresponding 
quantum map based on correlations between input and 
out states of quantum systems that serve as probes of the 
channel. 

In principle, the action of a quantum channel can be 
probed in two different ways: (1) The first option is to 
use as an input a single entangled state of a bi-partite sys- 
tem [1-4]. One particle (e.g., a qubit) of this bi-partite 
system is transformed under the action of the channel 
while the second particle remains unchanged (or evolves 
according to a known transformation). By performing a 
complete quantum tomography of the bi-partite system 
at the output of the channel and comparing the input 
and output states one can determine what is a specific ac- 
tion of the channel under consideration. (2) The second 
option is to use a collection of linearly independent test 
states (forming a basis of the vector space of all hermitian 
operators) [5-8]. By performing a correlation measure- 
ment between a specific input and corresponding output 
state (that has been tomographically reconstructed) we 
can determine the map that characterize the quantum 
channel. For a total determination of the map we have 
to use a complete set of test states. 

Reconstruction of quantum channels using entangled 
states might seem to be more efficient since a preparation 
of just one state of an entangled pair is required. Never- 
theless, there are two technical problems that make this 
approach less practical than a utilization of single particle 
states. Specifically, one has to generate input bi-partite 
entangled states with a very high fidelity (it is essential 
that for a reliable channel reconstruction the input test 
states have to be prepared with a very high fidelity). Si- 
multaneously, the tomographic reconstruction of a bipar- 



tite entangled state at the output of the channel has to 
be almost perfect. But the most difficult obstacle is to se- 
cure that the "reference" particle from the entangled pair 
does not undergo uncontrolled changes during the time 
when the second particle is affected by the action of the 
quantum channel. These conditions are rather difficult to 
met. For this reason in the present paper we will concen- 
trate our attention on the second scenario. This approach 
is based on the fact that the action of any channel is de- 
scribed by a linear map £ and therefore it is completely 
determined by its action on basis elements, i.e. a set of 
linearly independent states, which play the role of test 
states. We assume, that these single-partite test states 
are known (i.e., their preparation is under a complete 
control). Thus the process of the channel reconstruction 
reduces to the reconstructions of single-particle states at 
the output of the quantum channel. The number of test 
states equals dp , where d is the dimension of the Hilbert 
space of a quantum system under consideration. In order 
to fully characterize the action of the quantum channel 
acting on a such quantum system we need (P{d? — 1) real 
parameters, i.e. in the case of a qubit channel we need 
12 real parameters. 

In our previous paper [9] we have analyzed the question 
of the process reconstruction from incomplete experimen- 
tal data. That is, we have considered a situation when 
the correlation input-output measurements do not allow 
for a unique determination of the <P{(P — 1) real param- 
eters. We have shown how to proceed with the process 
estimation in the case of incomplete experimental data. 
In this paper we want to study more delicate problem - 
how to perform process estimation when the experimen- 
tal data are incomplete and/or they are not consistent. 
That is, the straightforward estimation leads to maps 
that are not physical (i.e., they are positive but not com- 
pletely positive) [10]. We will present a regularization 
procedure that allows us to handle such situations. More- 
over, we will show that this procedure also allow us to 
determine optimal approximations of non-physical oper- 
ations. As an example we will analyze in detail the so 
called universal NOT gate. The paper is organized as fol- 
lows. In Section 2 we briefly describe properties of qubit 
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state space while in Section 3 we discuss the structure 
of space of qubit operations (channels). Sections 4 and 
5 are devoted to the reconstruction of quantum channels 
from incomplete data. Approximations of non-physical 
operations are derived in Section 6. The conclusions are 
presented in Section 7. 



II. STRUCTURE OF QUBIT STATE SPACE 

Firstly, let us consider a simple geometrical represen- 
tation of a qubit state space. This state space has a 
topology of a sphere that is often called as the Bloch 
sphere. The set of all hermitian operators form a real vec- 
tor space endowed with a scalar product defined by the 
relation {A\B) = TrAB. Consequently, any operator can 
be written as a linear combination of operators that form 
an orthogonal basis. The set of the Pauli cr-matrices, i.e. 
{I , (Tx: <7y o'z} = {I^'^}: represents a standard choice of 
the operator basis for a qubit (sec, e.g., Ref. [8]). The 
state space of a qubit is a subset of all hermitian operators 
with a unit trace. Except the operator ao = I all other 
members of the <T-basis are traceless. Therefore any op- 
erator with a unit trace can be written as ^(Z + r-ir). 
Such operators represent a quantum state only when the 
operator q is positive, i.e. \f\ < 1. In this way we obtain 
the Bloch sphere representation of qubit states, g ^ f. 
A state reconstruction is a task of experimental specifica- 
tion of the vector f = {x,y, z). ^From the orthogonality 
condition TrcTfecr; = 2Ski we find that components of the 
vector r are determined by an expression fk = Trgak- 
That is, they are equal to mean values of the hermitian 
operators (measurements) (Tx,(7y,(7z- So the complete re- 
construction is straightforward: All one has to do is to 
measure mean values of three system operators cr,,. , ay , cr^ . 

Let us note that sometimes the reconstructed density 
operator may not satisfy the condition |r| < 1. This fail- 
ure of the reconstruction scheme is usually caused by an 
inconsistent measurement statistics which results in an 
incorrect identification of probabilities and consequently 
in the derivation of false mean values. The easiest way 
how to perform a regularization of the reconstruction in 
this case is the following one: The reconstructed density 
operator has to have always a unit trace, i.e. the operator 
is represented by a vector r though it might have a length 
larger than unity. In this case the reconstruction proce- 
dure fails since the estimated operator is not physical. 
One can argue that an actual physical state is the clos- 
est one to the "reconstructed" operator represented by a 
point on the Bloch sphere (a pure state) with rc point- 
ing in the same direction as the reconstructed vector r. 
Formally the regularization corresponds to a multiplica- 
tion of the original vector r by a positive constant k, i.e. 
Tc = kr. From the physical point of view this regular- 
ization can be understood as an admixture of a "white" 
noise described by the operator i/ [that is represented 
by the center of the Bloch sphere, i.e. = (0, 0, 0)] to the 



measured data. Formally this "regularization" procedure 

reads 

g, = kQ+{l-k)^I = ^{I + kf-a). (1) 

Such correction corresponds to the addition of completely 
random and equally distributed events ("clicks") to the 
outcome statistics of measurement results. In what fol- 
lows we will utilize an analogue of this intuitive picture 
to regularize reconstructions of maps describing quantum 
channels. 



III. STRUCTURE OF QUBIT CHANNELS 



The structure of qubit channels is known mainly due 
to work of M.B.Ruskai et al. [12]. Let us briefly summa- 
rize main properties of qubit channels. Any completely 
positive map £ can be imagined as an affine transfor- 
mation of the vector r, i.e. f r' = Tr + t. where 
T is a real 3x3 matrix and t is a translation. However, 
this form guarantees only the preservation of a trace and 
the hermiticity of the transformation £. In fact, the set 
of all completely positive tracepreserving maps forms a 
specific convex subset of all affine transformations. For 
qubits the number of parameters specifying the channel 
equals to 12. Because of the affinity of any evolution map 
£, one can use the following matrix representation 

The coefficients of the matrix £ arc given by a relation 
£ki = Tr(a/c£[c7i]), where (7fe(;) are Pauli a matrices. 

Any matrix T can be written in the so-called singular 

value decomposition, i.e. T = RijDRv where Ru,Rv 
are orthogonal rotations and D = diag{Ai, A2, A3} is a 
diagonal matrix with Afe being the singular values of T. 
Each three-dimensional orthogonal rotation Rjj (clement 
of the group S0(3)) is related to some qubit unitary trans- 
formation U (an clement of the group SU(2)) via the rela- 
tion U gU'^ = i(/ + {Rur) ■ a). This means that any map 
f is a member of less-parametric family of maps of the 
"diagonal form" <!>£. In particular, £[g] = U^£[V^ gYjU"^ 
where U, V are unitary operators. This reduction of pa- 
rameters is very helpful, and most of the properties (also 
complete positivity) of £ is reflected by the properties of 
$f . The map £ is completely positive (CP) only if 
is also CP. Let us note that is determined not only 
by the matrix D, but also by a new translation vector 
T = Rut, i.e. under the action of the map the Bloch 
sphere transforms as follows rj — > r'j = XjVj + tj. 
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FIG. 1. The set of positive unital trace-preserving linear 
maps "I>£ parametrized by three real parameters Ai,A2,A3. 
The point Ai = A2 = A3 = — 1 corresponds to a physi- 
cally unrealizable transformation - the universal NOT (de- 
noted as NOT in the figure) operation , i.e. The point 
Ai = A2 = A3 = 1 represents the identity map X (denoted 
as Id) 
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FIG. 2. The pictorial representation of the "position" of 
completely positive unital maps in the space of all positive 
unital maps. The CP unital maps form a tetrahedron with 
four unitary transformations in its vertices (extremal points) , 
that correspond to a-matrices. The un-physical universal 
NOT operation (denoted as NOT) and its best completely 
positive approximation called the optimal universal NOT (de- 
noted as U-NOT) are shown as well. The optimal universal 
NOT is represented by the point Ai = A2 = A3 = —1/3. 



A special type of completely positive maps is a class of 
unital transformations, i.e. maps for which the total mix- 
ture (the center of the Bloch sphere) is not affected by 
the transformation. It means that the translation term 
vanishes, i.e t = t = 0. In this case the geometrical 
analysis is quite simple. The positivity of the transfor- 
mation corresponds to the conditions |Afc| < 1, i.e. 
these transformations are represented by points lying in- 



side a cube (see Fig.l). The conditions of the complete 
positivity [12] requires the validity of the following four 
inequalities 

|Ai±A2| < |1±A3|. (3) 

These inequalities specify a tetrahedron inside a cube of 
all positive unital maps with the extreme points being 
four unitary transformations /, ax, c^, cr^ (see Fig. 2). As 
a result of this analysis one can conclude that the unital 
completely positive maps form a tetrahedron. 



IV. THE CHANNEL RECONSTRUCTION: 
INCONSISTENT STATISTICS 

Let us consider a situation when we want to estimate 
a channel, but neither experimental data nor non-trivial 
prior knowledge are available. Following the approach 
utilized for quantum state estimation [15] one can as- 
sume that in the absence of knowledge about the char- 
acter of quantum channel the most reliable estimation of 
this channel corresponds to an equally weighted "aver- 
age" over all possible quantum maps. The question is 
what is the average over all completely positive maps? 
We have already argued [9] that for qubits this average 
is the map A that transforms the whole state space into 
the total mixture, i.e. A[g] — \l ■ The reasoning goes as 
follows: The maps £± {T, ±t) have the property that 
if one of them is completely positive the second one is CP 
as well. Consequently, the average i {S+ + f - ) is a unital 
map, i.e. the average of all maps will be a unital map. As 
we already said, the unital maps (up to unitary transfor- 
mations) form a tetrahedron. The average over all points 
in tetrahedron is represented by the center of the tetra- 
hedron, i.e. Ai = A2 = A3 = 0, which corresponds to 
the contraction into the total mixture. Therefore, the 
average over all completely positive maps of acting on a 
qubit is the contraction of the whole Bloch sphere into 
its center, i.e. A[g] — \l. 

The result of a complete reconstruction based on the 
four test states (i.e., any collection of four mutually lin- 
early independent qubit states) is a map that is for sure 
trace-preserving and positive. However, it could happen 
that it is not completely positive. How to extract the 
physical map £c from an unphysical result £1 One way 
is to follow a similar reasoning like for the state recon- 
struction, when an unphysical result was corrected (reg- 
ularized) by adding a noise into the system. 

When the reconstructed map is not completely posi- 
tive we can regularize this result by an analogue of the 
total mixture, i.e. the map A. In particular, this regu- 
larization of a quantum channel reads 

^^ = k^Hl-k)A^{l/^^)- (4) 
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The correction (regularization) corresponds to a "mini- 
mal" adjustment of the parameter k such that the map 
£c is completely positive. 

Let us estimate what is this minimal (critical) value 
of the parameter k, i.e. the value which surely regular- 
izes any positive map and transforms it into a CP map. 
Trivially, it is enough to set fc = 0. In this case we 
completely ignore the measured data and the "corrected" 
map is A. However, wc are interested in some nontrivial 
bound, i.e. in the largest possible value of k. As we have 
already mentioned, the reconstructed map is always pos- 
itive. Consider, for simplicity, that the map is also uni- 
tal. Then the "worst" example of a positive map, which 
is not completely positive, is the universal NOT opera- 
tion. In this case the distance between this map and the 
tetrahedron of completely positive maps is extremal (see 
Fig. 2). This artificial example serves as a good test of 
our method, and gives us some bound on k, i.e. a value 
that surely corrects each result. The conditions of the 
complete positivity given in Eq.(3) imply that k = 1/3, 
i.e. Ai = A2 = A3 = —1/3 (see Fig. 2). Surprisingly, 
this is the same result as the one that has been obtained 
in Ref. [13] where the best (optimal) completely positive 
approximation of the universal NOT operation, i.e. an 
optimal universal NOT machine has been presented. In 
this sense our correction method works optimally. 



Becaiisc of the ad hoc assumption about the transfor- 
mation of unused test states it might happen that the 
resulting map is not completely positive. In this case 
the estimation procedure has to be complemented by a 
search for a map for which the total mixture is "shifted" 
as little as possible, i.e. the estimated channel is prefer- 
ably unital. A specific situation occurs when the data 
contain information about the transformation of the to- 
tal mixture. In this case the strategy suggests to use the 
state as the state on which all other undetermined 

states are mapped. As before, if such map is not com- 
pletely positive we have to search for a map, for which 
the deviation from the transformation £"[5/] is minimal. 

The method described above is discussed in detail in 
Ref. [9]. We note that this (incomplete) reconstruction 
can also fail (it gives no result), because it can happen 
that no physical channel is is compatible with given ex- 
perimental data. Even if the reconstructed operators g'j 
describe valid quantum states, the incomplete data can 
be in a contradiction with the; condition of a complete 
positivity. In what follows we will briefly describe our 
strategy on a particular example - estimation of the iden- 
tity channel (i.e. the channel, that does not change in- 
put states at all). For more detailed and more general 
description of this strategy see Ref. [9] . 



V. THE CHANNEL RECONSTRUCTION: 
INCOMPLETE DATA 

In this section we will present the reconstruction 
scheme which can be used when the number of test states 
is reduced so that the complete reconstruction of the 
channel cannot be performed [9]. On the other hand 
we assume that each of the state used is represented by 
an infinite ensemble of identically prepared states so the 
complete tomography of a corresponding state at the out- 
put of the channel can be performed. 

In the case of qubit channels the aim is to perform the 
reconstruction based on n = 0, 1,2, 3, 4 measured input- 
output correlations of the form Qj g'^ We have al- 
ready shown that having no information (n = 0), the 
best estimation of the map is the contraction into a to- 
tal mixture (an average over all quantum channels) i.e. 
£0 = A. Motivated by such result, the reconstruction 
strategy (sec Ref. [9]) is as follows: All undetermined 
states (belonging to the complement of the linear span 
of the used input states) are assumed to be transformed 
into the total mixture. That is, if a given test state is 
not explicitly used for a channel reconstruction (i.e., it 
is not know how this state is transformed by the action 
of the channel) it is assumed that the channel transforms 
this state into a total mixture. This additional assump- 
tion complements the knowledge of how other test states 
are actually transformed by the channel and allows us to 
use the deterministic procedure of channel reconstruction 
(for details see e.g. Refs. [6,5,7]). 



A. Case study: identity channel 

In what follows we will perform a step-by-step recon- 
struction of a qubit channel based on a knowledge of how 
a single, two and three test states are transformed under 
the action of a given channel. 

Single test state. Let us assume that our knowledge 
about the action of a particular channel is represented by 
the assignment 

Qi q'i = 81 ■ (5) 

We remind us that each state can be written in the 
form Qi = ^{I + f ■ a) = ^{I + wSz), where Sz = 
— |V'-l)(V'-l| with being eigenvectors of 

the operator q\ and w describes a n impurity of the state 
under consideration, i.e. w = -^/l — 2Tiq\. If w = 
then Qi describes maximally mixed state and for ii; = 1 
the state is pure. In the Bloch-sphere picture the pa- 
rameter w corresponds to the distance between the total 
mixture (the center of the sphere) and a point corre- 
sponding to the given state. One can define a new oper- 
ator basis Sx,Sy, Sz such that Sz is given as before, and 
Sj = UajU^ with unitary U. In this new basis the action 
of the channel is described by the matrix 

/ 1 \ 

c I X a d 

^'=\y b e 

\z c f ij 
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and the task of the estimation is to specify all matrix ele- 
ments. Our strategy suggests that all states (belonging to 
the complement of the linear span of gi ) are transformed 
into the total mixture, i.e. we set all the parameters equal 
to zero. 




FIG. 3. The pictorial description of the channel estima- 
tion when two states pi and p2 are used to test the action 
of a quantum channel and they are transformed according to 
Eq. (7). 

Two test states. In this case the knowledge about 
the action of the channel is represented by a transforma- 
tion of two test state given by the following relations 



Q2^ Q2- 



(7) 



Irrespective of states qi and Q2 we use, they specify a line 
crossing the Bloch sphere. In particular, from Eq. (7) we 
obtain a knowledge about the transformation of all states 
of the form 



g\ = >^gi + (1 - A)£i2 



(8) 



with a real parameter A. This line crosses the Bloch 
sphere in two extremal points representing pure states. 
Let us choose one of these pure states [denoted by = 
^{I+wSz)]- We also define a new operator (not necessar- 
ily positive) ^2 = ^{I + vSx), which in the Bloch-sphere 
picture corresponds to the intersection of the line g\ and 
the line orthogonal at the center of the Bloch sphere to 
the line given by points {£,1,^1} (see Fig. 3). This new 
state ^2 defines the operator Sx- Using the new operator 
basis Sx,Sy,Sz we can express the map characterizing 
the action of the channel as 



£2 = 



10 
X 1 a 
2/060 
z c 1 



(9) 



We again set all free (unspecified by the measurement) 
parameters to zero. As a result we find that the trans- 
formation £2 = diag{l, 1, 0, 1} is not completely positive, 
i.e. it does not describe a valid quantum channel. To pre- 
serve the unitality of the channel (f2[^-^] = ^I) one has 



to verify the complete positivity of the transformation 
£2 = diag{l, l,k,l}. We find out that the only possibil- 
ity is to take k = 1. Consequently, the identity channel 
is correctly estimated, but we are still not sure, whether 
the map is unital. However, we have to note one specific 
case, when the total mixture is among the states gx, i.e. 
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= i/. In this case 



/ 1 

e a 
0/50 

V 5 c 1 



(10) 



and our reconstruction procedure leads us to the channel 
£2 — diag{l, 0, 0, 1}, i.e. the values of all free parameters 
are equal to zero. Note that in this case £2 — £i- 

Three test states. Let us consider that the linear 
span of gi, g2, gz does not contain the total mixture. In 
this case the estimated transformation takes the form 



/ 1 

X 1 

y 1 

V z 1 



(11) 



The only possibility to preserve the complete positivity 
of £^ is to chose a; = y = z = 0, which is completely 
compatible with our strategy. Consequently, the channel 
is estimated perfectly and it is described by the trans- 
formation ^3 = diag{l, 1, 1, 1} = I. In the case, when 
from the measured data it follows that i/ i/, the es- 
timation coincides with the reconstruction with two test 
states, where the transformation of the total mixture is 
estimated to be £[\l] = \l- The only difference is that 
in the three-state case the unitality is guaranteed by the 
data. The reconstruction gives us the same result in both 
cases, i.e. £"3 =2. As a result we find that identity chan- 
nels for qubits can be uniquely identified using just three 
test states. 

Let us summarize the incomplete reconstruction of the 
identity channel. The hierarchy of estimations on dif- 
ferent levels specified by the number of test states is as 
follows 



£o^A; 

fi =diag{l, 0,0,1}: 
£2 = £3 — £4 — ^ ■ 



(12) 
(13) 
(14) 



/,From our previous discussion we can conclude that the 
identity channel can be reconstructed using just three 
test states. Given the fact, that for any unitary channel 
£ the induced map represents the identity channel, 
i.e. — 2 we can conclude that for a complete de- 
termination of a unitary channel we need just three test 
states. A unitary transformation is determined by the 
choice of the basis Sa:,Sy,Sz- 
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VI. COMBINATION OF DATA 
IMPERFECTIONS: INCOMPLETE AND 
INCONSISTENT DATA 

In this section we will study how to perform a process 
reconstruction when the data obtained from the measure- 
ment arc incomplete as well as inconsistent. To secure 
that the incomplete reconstruction does not fail we can 
adopt the regularization procedure as described earlier in 
the paper. In particular, we have two options: (i) either 
to regularize the output test states, or (ii) to regularize 
the estimated map itself. 

The first scenario has to be used always when the 
state reconstruction of some of the output test states 
fails. However, also if all the test states are esti- 
mated (reconstructed) correctly (i.e. they are k^gitimate 
physical states), the complete positivity is not guaran- 
teed and the process reconstruction might fail. This sit- 
uation can occur even for two test states. Specifically, 
the assignments q\ — > and qi — > ^2 ^"^^ compati- 
ble with some completely positive map E if and only if 
D{gi,tg2) > D{g[,tg'2) for all real positive t [14]. Unfor- 
tunately, for three and four test states no similar result is 
known. Therefore the complete positivity of the process 
reconstruction has to be checked for each individual case 
separately. 

In what follows we will utilize (demonstrate) two re- 
construction strategies: 

• We will use the reconstruction from incomplete data 

as discussed above. In the case of a failure (the 
reconstructed channel is not a CP map) a regular- 
ization of outputs of test states is performed (noise 
is admixed into the output states, so that the re- 
constructed map becomes CP). 

• We will assume that an unknown map transforms all 

states except the basis test states into the total mix- 
ture, i.e. we set all free parameters to zero. If the 
resulting map is not completely positive, then the 
map will be regularized by adding the average chan- 
nel A. 

Let us consider a particular example that allows us 

to demonstrate these two methods. In our example 
we will use an "artificial" and unphysical data gener- 
ated by the universal NOT operation [13]. The logi- 
cal NOT operation is defined by relations |0) — » |1), 
|1) |0) in a computer basis {|0);|1)}. These rela- 
tions do not completely determine a quantum channel, 
i.e. many quantum channels perform such transforma- 
tions. For instance, the "classical" NOT can be viewed 
as the transformation of the Bloch sphere into the line 
connecting North (|1)(1|) and South Pole (|0)(0|), i.e. 
mTe[g] = ^{I + {(Jz)eS,) with {cj,)g = Tvga,. There 
exists also a unitary ("quantum") realization of the log- 
ical NOT operation, i.e. NOTg[^p] = Uxgc^x- However, a 
natural generalization of the NOT operation is the so 



called universal NOT gate, which performs the transfor- 
mation 1-0) for all states j-f/;). The universal NOT 
gate is not completely positive (for mixed states we will 
use the notation g g'^). The optimal universal NOT 
operation is the closest physically valid map that per- 
forms approximatively the universal NOT. In Ref. [13] it 
has been shown that the corresponding quantum chan- 
nel that maximizes the average fidelity under given con- 
straints reads ^not = diag{l, —1/3, —1/3, —1/3}. 

A. Case study: the universal NOT gate 

Formally the universal NOT operation determines the 
following transformations on a set of test (basis) states 



Qx = 




^ Q'x = 


\ii-(^x) ; 


(15) 


Qy = 


\{I + ^y) 


^ q'v = 


\ii-'^y) ; 


(16) 


Qz = 








(17) 




Qo = \l 






(18) 



Let us start with the first reconstruction strategy as de- 
scribed above. In the case when just one or two states 
have been used to test the action of the quantum channel 
no regularization is needed to estimate the channel. 

(1) When we use only a single test state (e.g., a pmc 
state) the resulting estimated map is a contraction of the 
Bloch sphere into a line connecting two mutually orthog- 
onal states. If we use the data g^ —>■ gl^, then this line 
is given by the points (0, 0, ±1), i.e. the Bloch sphere is 
mapped into the z axis. This map can be understand as 
the "classical" logical NOT. 

(2) In the case of two test states, we obtain a unitary 

rotation. If we use g^ (J^ ^-nd gy g'y, then the result 
is a rotation by the angle tt around the x axis, i.e. the 
(Tx operation. This operation is usually referred to as the 
"quantum" (unitary) logical NOT (see Fig. 4). 

When three test states are used to determine the ac- 
tion of the universal NOT gate we face a serious diffi- 
culty: Let us consider transformations of three test states 
Qx,y,z Q'x.y,z- t^is casc uo Completely positive map 
exists. Therefore, a regularization of output states is 
required - in particular, we have to determine a mini- 
mal amount of noise that is "included" in the output 
states so that the channel estimation with these regular- 
ized (noisy) output states will result in a CP map. It 
turns out that the amount of noise corresponds to the 
value k = 1/3. Therefore the transformation that is clos- 
est to the universal NOT gate and is CP has the form 



We can conclude that already for three test states we 
find that the reconstructed map that is CP is the opti- 
mal universal NOT operation [13]. In Fig.4 we see how 
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the estimation of the map changes with the number n of 
used test states. 




FIG. 4. The figure represents how our knowledge about the 
channel (the universal NOT gate) is improved with increasing 
the number of test states when the first reconstruction strat- 
egy is used. If no measurement has been performed then the 
best possible estimation of the channel is the contraction of 
the Bloch sphere into its center (n = 0). When a single test 
state Qz q'z is used the optimal estimation of the trans- 
formation is a contraction of the whole Bloch sphere into a 
line connecting the "north" (denoted as A*') and "south" poles 
of the Bloch sphere (i.e. the points (0,0, ±1), respectively). 
This contraction is represented by the figure n = 1 in the 
picture. When two test states (n = 2) are used, i.e. Qz q'z 
and Qy — > Qy, then the result of the estimation of the action 
of the quantum channel is a rotation by the angle vr around 
the X axis, i.e. the operation. Finally, when three (n — 3) 
or four (n = 4) test states are used to determine the action of 
the quantum channel prescribed by the universal NOT trans- 
formation, the best estimation of the channel is the optimal 
universal NOT gate (see the two figures of the lower line in 
the picture, both figures are the same and they describe the 
action of the optimal universal NOT gate). In this case the 
regularization of the output test states has been used in order 
to estimate a map that is completely positive. As a specific 
example we show how the "north" pole N of the Bloch sphere 
is transformed under the action of estimated maps. 



In what follows we will study the second strategy for a 
reconstruction (estimation) of the best possible approxi- 
mation of the universal NOT gate. Instead of searching 
for a completely positive map by "admixing" the small- 
est possible amount of noise into the outputs of the test 
states, let us assume that all states except the test states 
are mapped into the total mixture ^ . Using this approach 
we can be sure that the resulting map is positive, but not 
completely positive. Therefore we have to use a regular- 
ization of the map. Obviously, in this case the resulting 



map will not explicitly satisfy conditions imposed by the 
transformation of test states, Qj — > g'j, i.e. with the given 
data. But this has to be expected since it is the only way 
how to impose the CP condition on the reconstructed 
map. Using this type of regularization of the map we 
arrive at the same estimations as in the case of previous 
strategy. However, for n = 2 (two test states) the situa- 
tion is different. In this case, the reconstructed map acts 
as follows 

f=ix,y,z)^f' = iO,-y,-z). (20) 

That is, the states of the form g — ■^{I + xcjx) are trans- 
formed into the total mixture. Consequently, the map is 
unital and has a diagonal form with Ai = 0, A2 = A3 = 
— 1, which is not compatible with the complete positivity 
(see Fig. 2). The parameter fc = 1/2 can be used to cor- 
rect this map. This value of k can be derived from the 
inequalities |Ai ± A2I < |1 ± A3I. One can observe this re- 
sult in Fig. 2 as the closest point from the tetrahedron to 
the point (0, —1, —1). The whole incomplete reconstruc- 
tion using this method is depicted in Fig. 5. In Tab. 1 
we present diagonal elements (all other matrix elements 
are equal to zero) of matrices corresponding to recon- 
structed maps based on the results of measurement of n 
test states. We present results for both reconstruction 
methods. 





1st method 


2nd method 




{1,0,0,0} 


{1,0,0,0} 


£1 


{1,0,0,1} 


{1,0,0,1} 


£2 


{1,1,-1,-1} 


{1,0,-1/2,-1/2} 


£3 


{1,-1/3,-1/3,-1/3} 


{1,-1/3,-1/3,-1/3} 


£4 


{1,-1/3,-1/3,-1/3} 


{1,-1/3,-1/3,-1/3} 



TABLE I. The two reconstruction scenarios for an 
un-physical universal NOT gate result in a sequence of maps 
depending on the number of test states. The subscript in the 
description of a given map £j corresponds to the number of 
test states that have been used in the reconstruction. From 
the table it is clear, that both methods give us the same pro- 
cess estimation except for the case when two test states have 
been used. 



^More specifically, states that are complemented to a linear span of the test states. 
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FIG. 5. The figure represents how our knowledge about the 
channel (the universal NOT gate) is improved with increas- 
ing the number of test states when the second reconstruction 
strategy is used. If no measurement has been performed then 
the best possible estimation of the channel is the contrac- 
tion of the Bloch sphere into its center (n — 0). When a 
single test state Qz —* g'^ is used the optimal estimation of 
the transformation is a contraction of the whole Bloch sphere 
into a line connecting the "north" and "south" poles of the 
Bloch sphere (i.e. the points (0,0, ±1), respectively). This 
contraction is represented by the figure n = 1 in the picture. 
When two test states (n = 2) are used, i.e. — > g'^ and 
Qy ~* Qy^ then the result of the estimation of the action of 
the quantum channel is described by the transformation (19). 
Consequently, the map is unital and has the diagonal form 
with Ai = 0, A2 = As = —1. Unfortunately this map is not 
completely positive. The regularization parameter k — 1/2 
is used to correct this map. Finally, when three (n — 3) or 
four (n = 4) test states are used to determine the action of 
the quantum channel prescribed by the universal NOT trans- 
formation, the best estimation of the channel is the optimal 
universal NOT gate. 



can be understood as an addition of a white noise into our 
data. In principle, we can face two situations: Either re- 
constructions of all test states at the output of the chan- 
nel correspond to proper quantum states, or there exist 
some outcomes that are not proper quantum states. In 
this case each of these states can be corrected by adding 
some noise (using the multiplicative factor k as discussed 
above). We have to keep in mind that the parameter k 
can be state dependent. Therefore, we have to choose 
the smallest value of k in order to correct the map. Once 
all the outputs (as estimated from the measured data) 
are proper physical states we can start to reconstruct the 
map itself. In spite of the fact that all the test states at 
the output are proper physical states the map that is esti- 
mated on the basis of these states may be not completely 
positive. In this case we have to search for the largest k, 
for which the corrected map £c = k£ + {1 — k)A is com- 
pletely positive. We have shown that at least for unital 
maps the value k — 1/3 always regularizes the estimated 
map. 

In order to illustrate our methods we have used the 
data generated by an unphysical map - the universal 
NOT gate. We have shown that using both approaches 
the result is the same and the reconstructed map is the 
best (optimal) approximation of the NOT operation [13]. 
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VII. CONCLUSION 



Reconstruction of quantum maps is a challenging prob- 
lem motivated mainly by experimental realizations of 
quantum gates. These gates have to be tested thoroughly 
in order to use them for any quantum computation. An- 
other important application of a channel reconstruction 
is in quantum communication when characteristics of a 
quantum channel have to be determined from a limited 
set of tests performed on the channel. 

In the present paper we have discussed some strate- 
gies how to estimate quantum channels when only in- 
complete and/or incompatible data from measurements 
is available. The incompatibility of available data results 
in estimated maps that are not completely positive. In 
this case a regularization of the reconstruction procedure 
is required in order to recover a physical (CP) map. 

The regularization method represents a correction of 
insufficient statistics. It uses a single parameter k that 



[1] G.M. D'Ariano and P.Lo Presti, Tomography of quantum 
operations, Phys.Rev.Lett. 86, 4195 (2001). 

[2] G.M. D'Ariano and P.Lo Presti, "Characterization of 
quantum devices" in Quantum Estimations: Theory and 
Experiment, Springer Series on Lecture Notes in Physics, 
Band 649, eds. G.M.Paris and J. Rehacek (Springer- 
Verlag, Berlin, 2004), p. 299. 

[3] M.Jezek, J.Fiurasek, Z.Hradil, "Quantum inference of 
states and processes", Phys.Rev.A 68, 012305 (2003). 

[4] Z. Hradil, J. Rehacek, J. Fiurasek, and M. Jezek, 
Maximum-Likelihood Methods in Quantum Mechanics in 
Quantum State Estimation, Springer Series on Lecture 
Notes in Physics, Band 649, eds. M.G.A. Paris and 
J. Rehacek (Springer- Verlag, Berlin, 2004), p. 63. 

[5] J.F. Poyatos, J.I. Cirac, and P. ZoUer, Complete charac- 
terization of a quantum process: two-qubit quantum gate, 
Phys. Rev. Lett. 78, 390 (1997). 



8 



[6] I.L. Chuang and M.A. Nielsen, Prescription for crperi- 

mental determination of the dynamics of a quantum black 

box, J. Mod. Opt. 44, 2455 (1997). 
[7] V. Buzek, Reconstruction of Liouvillian superoperators, 

Phys. Rev. A 58, 1753 (1998). 
[8] M.A. Nielsen and I.L. Chuang, Quantum Computation 

and Quantum Information (University Press, Cambridge, 

2000). 

[9] M. Zinian, M. Plescli, and V. Buzek, On the reconstruc- 
tion of superoperators from incomplete data, to appear in 
Foundations of Physics, qu£ait-ph/0406088. 

[10] M. Ziman, T. Hannemann, M. Plesch, V. Buzek, and Ch. 
Wunderlich, in preparation 

[11] M. Ziman, M. Plesch, and V. Buzek, in preparation 



[12] M.B. Ruskai, S. Szarck, and E. Werner, An analysis of 
completely positive tracepreserving maps on 2x2 matrices, 
Lin. Alg. Appl. 347, 159 (2002). 

[13] V. Buzek, M. Hillery, and R.F. Werner, Optimal manip- 
ulation with qubits:Universal-NOT gate, Phys. Rev. A 60, 
R2626 (1999). 

[14] A. Carlini and M. Sasaki, Geometrical conditions for 
completely positive trace-preserving maps and their appli- 
cation to a quantum repeater and a state- dependent quan- 
tum cloning machine, Phys. Rev. A 68, 042327 (2003). 

[15] V. Buzek and R. Dorka, Quantum Observations, in Co- 
herence and Statistics of Photons and Atoms, ed. J. 
Pefina (John Wiley & Sons, New York, 2001), p. 198. 



9 



